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ABSTRACT 



Dynamical properties of two-component galaxy models whose stellar density distribu- 
tion is described by a 7-model while the total density distribution has a pure r^^ profile, 
are presented. The orbital structure of the stellar component is described by Osipkov-Merritt 
anisotropy, while the dark matter halo is isotropic. After a description of minimum halo mod- 
els, the positivity of the phase-space density (the model consistency) is investigated, and nec- 
essary and sufficient conditions for consistency are obtained analytically as a function of the 
stellar inner density slope 7 and anisotropy radius. The explicit phase-space distribution func- 
tion is recovered for integer values of 7, and it is shown that while models with 7 > 4/17 
are consistent when the anisotropy radius is larger than a critical value (dependent on 7), the 
7 = models are unphysical even in the fully isotropic case. The Jeans equations for the stel- 
lar component are then solved analytically; in addition, the projected velocity dispersion at 
the center and at large radii are also obtained analytically for generic values of the anisotropy 
radius, and it is found that they are given by remarkably simple expressions. The presented 
models, even though highly idealized, can be useful as starting point for more advanced mod- 
eling of the mass distribution of elliptical galaxies in studies combining stellar dynamics and 
gravitational lensing. 
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1 INTRODUCTION 

Analysis of stellar kinematics (e.g. Bertin et al. 1994, Rix et al. 
1997, Gerhard et al. 2001), as well as several studies combining 
stellar dynamics and gravitational lensing strongly support the idea 
that the dark and the stellar matter in elliptical galaxies are dis- 
tributed so that their total mass profile is described by a density 
distribution proportional to (e.g., see Treu & Koopmans 2002, 
2004; Rusin et al. 2003; Rusin & Kochanek 2005; Koopmans et al. 
2006; Czoske et al. 2008; Dye et al. 2008). In particular, Gavazzi 
et al. (2007), with a gravitational lensing analysis of 22 early-type 
strong lens galaxies, reported a total r^^ density profile in the range 
1-100 effective radii. It is clear that in this field the availability of 
simple dynamical models of two-component galaxies can be use- 
ful as starting point of more sophisticated investigations based on 
axysimmetric or triaxial galaxy models (e.g., Cappellari et al. 2007, 
van den Bosch et al. 2008). A few simple yet interesting models 
with flat rotation curve have been in fact constructed, such as those 
in which the stellar mass was described by a power-law in a to- 
tal mass distribution (e.g. Kochaneck 1994), or those obtained 
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from physical arguments (in case of disk galaxies, see e.g. Naab & 
Ostriker 2007). 

Here the family of two-component galaxy models whose total 
mass density is proportional to r~^, while the visible (stellar) mass 
is described by the well-known 7 models (Dehnen 1993, Tremaine 
et al. 1994), is presented. Some preliminary numerical investiga- 
tion of these models has been done in Keeton (2001), and they 
have been used in Nipoti et al. (2008) as diagnostics of the total 
mass distribution in elliptical galaxies. In this paper a more system- 
atic study of the dynamical properties of these models is presented. 
It is shown that the Jeans equations for the stellar component with 
Osipkov-Merritt (Osipkov 1979, Merritt 1985, hereafter OM) radial 
anisotropy can be solved analytically. Remarkably, the projected 
velocity dispersion at the center and at large radii can be expressed 
in terms of the model circular velocity by means of extremely sim- 
ple formulae for generic values of the anisotropy radius and of the 
central stellar density slope 7. In principle this feature opens the 
possibility to obtain preliminary indications about the anisotropy 
from observations at small and large radii. The positivity of the 
phase-space density (the so-called consistency) is investigated, by 
obtaining analytically the necessary and sufficient conditions for 
model consistency in terms of 7, of the anisotropy radius, and of the 
dark- to-stellar mass ratio within some prescribed radius. It is found 
that the phase-space distribution function (hereafter DF) can be re- 



2 Ciotti, Morganti, & de Zeeuw 



covered analytically for 7 = 0, 1, and 2. In particular, it is shown 
that 7 = models in a total density profile are unphysical for 
any value of the anisotropy radius. These results extend the class of 
two-component galaxy models with explicit DF and add to the large 
amount of phase-space information already available about one and 
two-component 7 models (e.g., see Dehnen 1993, Tremaine et al. 
1994, Hiotelis 1994, CaroUo et al. 1995, Ciotti 1996, 1999; Baes et 
al. 2005, Buyle et al. 2007, Ciotti & Morganti 2008). 

The paper is organized as follows. In Section 2 the main struc- 
tural properties of the models are presented, while in Section 3 an 
investigation of the phase-space properties of the models is car- 
ried out both from the point of view of necessary and sufficient 
conditions for consistency and of direct recovery of the DF in spe- 
cific cases. In Section 4 the solution of the Jeans equation with 
OM radial anisotropy is presented, together with their projection at 
small and large radii. Finally a short summary of possible use of 
the present models in observational works is given. 



2 THE MODELS 

2.1 Stellar distribution 

The density profile of spherical 7 models is 

- (3 - 7)M, 



p*(r) = 



s~i{l + sY-~i ■ 



A, 



4:-Kri 



(1) 



where < 7 < 3, M, is the total stellar mass, r, is a scale-length, 
and s = r/r* is the dimensionless radius. These models have been 
investigated extensively, and here only the properties of present use 
are listed. In particular, the cumulative stellar mass within r is given 
by 



M,{r) = M. 



3-7 



(2) 



so that the dimensionless half-mass spatial radius is Sh ~ 
1/(23-"' — 1). The projected stellar surface density 

pt{r)rdr 



E.(i?) = 2 



(3) 



(e.g., Binney & Tremaine 2008) cannot be expressed in terms of el- 
ementary functions for generic values of 7, however the asymptotic 
behaviour of E*(i?) is easily obtained for i? — > 

^ (0<7<1); 

(7 = 1); (4) 
(1<7<3); 



(3-7)(2-7)(l-7)' 
-Slogry, 

0Fr(7/2 - 1/2) 



and for J? ^ 00 



2r;3 



r(7/2) 



(0 < 7 < 3). 



(5) 



In the equations above R is the radius on the projection plane, 
r] = i?/r« is its normalized value, and F is the complete Euler 
Gamma function; note that the first of eqs. lHJl is the exact value of 
the projection integral for i? = when < 7 < 1. 



2.2 Total and dark matter distribution 

By assumption the total mass density is taken to be 



PT{r) 



TIA, 



(6) 



where 7?. is a dimensionless scale factor which measures the im- 
portance of the dark matter density with respect to the stellar one: 
therefore, the stellar distribution would be a tracer in the total den- 
sity distribution in the formal limit A, — > and TZ — > 00, in a way 
such that the product TZ-A^ remains constant. The cumulative total 
mass within r is 



MT(r) = ATiTlA,rl s, 



(V) 



and the system (constant) circular velocity is — 47vGTZA,r:,; 
from this expression the dimensionless constant TZ (or the density 
scale j4*) everywhere it appears in favor of Vc- The total projected 
mass density at R is obtained from eqs. l|3) and l|6j as 



ET(i?) 



■KTZA,r, 



(8) 



so that the total mass contained within the cylinder of radius R is 

MpTiR) = 2-n I T.T{R)RdR = 27^^7^A.r2^?. (9) 
Jo 

Not all values of the coefficient TZ and of the inner stellar density 
slope 7 are compatible. In fact a first limitation is given by the 
request of positivity for the halo density 



, , A, 



7^- 
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This request restricts the value of 7 to the interval < 7 < 2, 
independently of the value of TZ. With 7 in the acceptable range, 
Ph is positive provided that 



7^ > 7^™,(7) = 



(4-7) 



4 — 7 



(11) 



(see Appendix A). For example, TZm{0) ~ 1/16, 7?.m(l/2) = 
0.0916, 7^m(l) = 4/27, and7^m(2) = 1; in Fig.[I](bottom panel) 
the minimum value TZm{'y) for halo positivity is represented by the 
solid line. A dark halo with TZm is called a minimum halo. While 
the density distribution of the minimum halo increases at the center 
as r^"^ for < 7 < 2, for 7 = 2 it results ph oc r^^, and so 
minimum halo 7 = 2 models are more and more baryon dominated 
near the center. We remarlt that the local mass-to-light ratio, 
proportional to p-rir) /p, (r) under the hypothesis of a constant 
stellar mass-to-light ratio, is a non-monotonic function of r as 
it increases near the center and for r ^ 00. The only exception 
is represented by the 7 = 2 case, which is characterized by a 
monotonically increasing mass-to-light ratio for increasing r. 

Of course, the positivity of ph is just a first condition for the 
acceptability of the model. A plausible second request is the mono- 
tonicity of ph as a function of radius; while at this stage monotonic- 
ity reduces to the determination of a minimum value of TZm (7) so 
that dph/dr < 0, in Section 3 it will be shown that this request is 
based on deeper physical arguments than simple structural plausi- 
bility. The explicit calculation of this additional restriction of TZ is 
given in Appendix A, and the resulting function TZm (7) is shown 
in Fig. [T] (bottom panel) with the dotted line: it is apparent that the 
request of monotonicity is just a little bit more stringent than posi- 
tivity, and that in the 7 = 2 case the two requests coincide. 

It can be of interest in applications to evaluate the relative 
amount of dark to visible mass within a prescribed radius. This 
quantity is easily calculated from eqs. ^ and Q. For example, 
within the half-mass radius rh one has 
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Figure 1. Bottom panel: minimum value of 7?. as a function of 7 for the halo 
density positivity (solid line, eq. 1111 ). monotonicity (dotted line, eq. lAll l 
and to satisfy the WSC in the isotropic case (dashed line, eqs. 1181- IA4I ). 
Top panel: the dark to stellar mass ratio within the half-mass radius of the 
stellar component, as given in eq. (12), for the three Umits in the botton 
panel. 



Mh(rh) ^ 2(3-7)7^™(7) 



21/(3-7) 



1, 



(12) 



where Mh (r) — A/t (r) — M, (r ) . In Fig.[T](top panel) the mass ra- 
tios corresponding to the three limits in the bottom panel are shown. 
The values are smaller than unity for all values of 7, except for the 
7 = 2 case. 

Another observationally relevant quantity is the projected 
mass ratio of dark-to-visible, for example within the effective ra- 
dius Rc of the stellar distribution. This is given by 



2Mph(i?e) 



> 7r(3 - ■y)'Jlm{j)'no 



Rc/r 



(13) 



where for example rjc ~ 2.9036 for 7 = (Dehnen 1993), rjc — 
2.3585 for 7 = 1/2, rye ^ 1.8153 for7 = 1 (Hernquist 1990), and 
r/c ~ 0.7447 for 7 = 2 (Carollo et al. 1995; note that in Jaffe 1983 
the slightly erroneous coefficient 0.763 is reported). These values 
translate into mass ratios of ~ 0.71, 0.70, 0.69, and 1.34 when 



considering the minimum value of TZ,n{'y) for halo positivity, for 
7 = 0, 1/2, 1, and 2, respectively. 



3 THE PHASE-SPACE DISTRIBUTION FUNCTION 

Before solving the Jeans equations, it is useful to discuss some ba- 
sic property of the phase-space DF of the presented models, in or- 
der to exclude physically inconsistent combinations of parameters 
(i.e., choices that would correspond to a somewhere negative DF). 
Fortunately, as discussed extensively in Ciotti & Pellegrini (1992, 
hereafter CP92), C96, and C99, it is possible to obtain lower bounds 
for the CM anisotropy radius as a function of the density slope and 
the total mass profile, without actually recovering the DF, which 
is in general impossible in terms of elementary functions. More 
specifically, in CP92 a simple theorem was proved regarding the 
necessary and sufficient limitations on ra in multi-component CM 
models, while more recently An & Evans (2006) proved the so- 
called "cusp slope-central anisotropy" theorem (see also eq. [28] in 
de Bruijne et al. 1996): the link between the two results is briefly 
addressed in Ciotti & Morganti (2008). Before using the CP92 test 
some preliminary work is however in order, because at variance 
with the common case of finite total mass, the total potential 



= Vr Ins 



(14) 



is now quite peculiar, being logarithmic. This means that in princi- 
ple orbits of any energy can be present, and the standard OM pre- 
scription must be reformulated to take into account the divergent 
behaviour of the potential both at r = and r = 00. Thus, a DF 
with the functional dependence 



f = f{Q), 



2ri 



(15) 



is assumed, where E — (f>T+v^ /2 and J are the energy and angular 
momentum modulus of each star (per unit mass), respectively. Note 
that, at variance with the usual OM parameterization, no cut on / 
for negative Q is present. By integration over the velocity space it 
is easy to show that for a given density component (stars or halo) 
in the total potential ^t, the density is related to its DF by 



47r 



1 + r^/r 



(16) 



</>T 



in principle, can be different for stars and dark matter. The anal- 
ogy with the standard OM relation is apparent (e.g., see Binney 
& Tremaine 2008). Following a similar treatment, it can also be 
shown that the radial ((Jr) and tangential (at) components of the 
velocity dispersion tensor are related as in the standard OM case, 
i.e. 

a?(r) 



2a? (r) 



(17) 



so that the fully isotropic case is obtained for 00, while for 

Ta = the galaxy is supported by pure radial anisotropy. For finite 
values of ra, the velocity dispersion tensor becomes isotropic for 
r (in practice for r < ra), and fully radially anisotropic for 
r ^ cx) (in practice for r > ra). Introducing the augmented density 

g{r) = p{r){l + i,], (18) 



eq. ( 116b can be Abel inverted, obtaining 

1 d f°° dg d(f>T 



/(0) = 



VStv^ dQ 7„ d(j}T 
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d? Q d4>T 



(19) 



where it is intended tiiat g is expressed in terms of (j>T, and tlie 
second identity follows from integration by parts when considering 
the untruncated nature of the studied density distributions. 

Moreover, for the present class of models it can be also proved 
that the velocity profile (VP, e.g. Carollo et al. 1995) can be written 



E.VP = 47r 
with 

g{r,R) 



g{r, R)rdr 



f{Q)dQ, 



^/r^ + r'iy/r^ + r'i 



i?2 



and 



= 0T + 



2 , 2 



'II 



7?2 2 ' 



(20) 



(21) 



(22) 



where U| | is the velocity along the line of sight direction. The inner 
integral in eq. l |20t can be simplified by using the first identity in 
eq. 



3.1 Necessary and sufficient conditions for consistency 

Repeating the same treatment of CP92, after differentiation of 
eq. ( I16t with respect to (j>T, it follows that a necessary condition 
for the positivity of the DF is that 



<0 [NC]. 



(23) 



de(r) 
dr 

This necessary condition for the DF positivity is independent of the 
radial dependence of the other density components of the system. 

The CP92 weak sufficient condition for consistency is recov- 
ered by requiring that the second derivative inside the integral in 
eq. M9\ be positive. In analogy with eq. J23b . this condition can be 
expressed as a function of radius as 



d_ 

dr 



dQ(r) 



> [WSC], 



(24) 



dr MT(r) 
where the total mass is given by eq. 

3.1.1 Isotropic halo consistency 

The first application of eqs. l l23t and l l24t concerns the consistency 
of the halo density distribution ph. For simplicity we restrict to the 
isotropic case, and then eq. ( I23t shows the equivalence of the re- 
quest of monotonicity of ph discussed in Sect. 2.2 with the nec- 
essary condition for a phase-space consistent halo. The WSC for 
a fully isotropic halo can be discussed analytically as described in 
Appendix A, and the resulting limit is represented by the dashed 
line in Fig.l (bottom panel): note how the three conditions of 
halo positivity, monotonicity, and consistency produce very simi- 
lar curves, that coincide for 7 = 2. 

Of course, tlie restriction of tfie study to an isotropic lialo is 
quite arbitrary, as tlie virialized end-states of A'^-body collapses 
are characterized by some amount of orbital anisotropy (e.g., 
van Albada 1982; Nipoti, Londrillo & Ciotti 2006). However, 
the present investigation is mainly focused on the properties of 
the visible component, and therefore we adopt the simplest dy- 
namical structure for the halo. 




Figure 2. The NC limit for consistency of 7 models is shown with the 
dotted line: all values of the anisotropy radius below the line con'espond 
to inconsistent models (C99). Solid line: the WSC for the present models, 
i.e., the locus above which the models are certainly consistent. The line 
begins at 7 = 4/17, as the WSC is not satisfied for centrally flatter models. 
Dashed line: the WSC for one-component 7 models as determined in C99. 
The triangles are the true anisotropy limits for the 7=1 and 7 = 2 models 
with halo obtained from their DF (Section 3.2). 



3. 1.2 Anisotropic stellar consistency 

In general, when dealing with OM anisotropic systems, the inves- 
tigation of the NC and WSC (eqs. (23l-(24l), and the study of the 
DF positivity (eq. 1191 ) lead to inequalities of the kind 



G 



> 0, 



(25) 



that must hold over all the domain of interest, which we indicate 
with C. In practice, the functions F and G are radial functions (in 
the case of the NC and WSC) or functions of the phase-space vari- 
able Q (in the case of the DF), depending on the specific context. 
Following C99, here we recall that according to eq. l |25| l all OM 
models can be divided in two main families. In the first case, the 
function F is nowhere negative over C (this could be the case of a 
system with a positive isotropic DF). Then, the lower hound to Sa 
is given by 



0,supc (j-^) 



(26) 



and the condition J25l > is satisfied provided Sa, > s^: in particular, 
if G is also positive over C then the system can be supported by 
radial orbits only. However, it may happen that F is positive only 
over some proper subset C+ of C, and negative (or zero) over C- . It 
trivially follows that in this second class of models if G is negative 
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over some subse{3 of C_, then the condition l l25b cannot be satis- 
fied for any value of Sa. On the contrary, it may happen that G is 
everywhere positive on C- : in this case one must consider not only 
the lower limit over C+, but also the upper bound 



infc 



(27) 



so that Sa < Sa over C_ for consistency. Summarizing, if F > 
over all its domain (i.e., C+ = C), then Sa > s~ satisfies inequal- 
ity ^5i. If F < over some set C_ but G > there, then the 
inequality s" < Sa < s+ must be verified. Finally, if over C- 
Sa < Sa or G < somewhere, then inequality l |25l l cannot be 
satisfied and, in case of a DF analysis, the model must be rejected 
as inconsistent. 

For example, in the case of 7 models the Sa (7) limit from 
the NC has been calculated analytically in C99, and the critical 
value of the anisotropy radius expressed in units of r* are ~ 0.354, 
~ 0.128, and for 7 = 0,1,2, respectively. In other words, smaller 
values of Sa correspond to physically inconsistent models (even 
though the solution of the associated Jeans equations is positive 
- see the following Section). The anisotropy limit over the whole 
range of 7 here considered is represented with the dotted line in 
Fig. 2: note how centrally flatter models are associated with larger 
values of the critical anisotropy radius. 

We now discuss the case of the WSC for the stellar component 
of our models. As shown in Appendix A, simple algebra reveals 
that the function F in equation l |25t is positive everywhere for 7 > 
4/17, and from eq. iA5\ the maximum of the function —G/F can 
be determined by solving an equation of degree four. The resulting 
value of Sa (7) is shown with the solid line in Fig. [2] For < 
7 < 4/17 the function F has two positive roots, delimiting the 
interval C- on which F < 0. The function G is positive on C- , so 
that we can determine the two values Sa (on C-) and Sa (on C+): 
however it turns out that Sa < Sa for < 7 < 4/17, so that the 
WSC is not satisfied, and for this reason the solid line interrupts in 
Fig.|2] Of course, being just a sufficient condition, this result does 
not exclude that consistent models exist for 7 < 4/17, but this is 
not assured as it is for the models with 7 > 4/17. For reference, 
in Fig.[2]the dashed line represents the WSC for one-component 7 
models as derived in C99. In particular, note how for models with 
7>1, the presence of the halo appears to increase the model ability 
to sustain radial orbital anisotropy, while flatter models in presence 
of the halo are less able to sustain anisotropy. 



3.2 Explicit Phase-Space DF 

For generic 7, eq. l |I9l ) can be rewritten as 

A, d f°° dg d<f 



A, 



(28) 



where '3/ = 4>i/v1 and q = Q/v^. The function g is the aug- 
mented density in eq. (18) normalized to A,, expressed in terms 
of the total potential. This is accomplished by elimination of the 
radius from the dimensionless identity s = exp(\I') obtained from 



^ In C99 and Ciotti (2000) it is eiToneously stated that the model is incon- 
sistent if G < everywhere on C_ . All the results presented therein are 
however correct. 




Figure 3. The phase-space DF (normalized to /v^\/?tTT'^) of the stellai' 
component of 7 = 1 (top) and 7 = 2 (bottom) models embedded in a dark 
matter halo so that the total density profile is proportional to . Solid 
lines refer to the case of a fully isotropic stellar component, dotted lines to 
intermediate values of the (nomialized) anisotropy radius Sa (1 for 7 = 1 
and 0.1 for 7 = 2), and finally the dashed lines to a value of Sa very near 
to the critical value for consistency. 



eq. (14). Not surprisingly, for generic 7 the functions U and V can- 
not be expressed in terms of known functions, however in Appendix 
B it is shown that for 7 = 0,1 and 2 the functions U and V can be 
expressed as simple linear combinations of exponentials and Poly- 
logarithms. Numerical inspection of the DFs shows that for 7 = 1 
and 2 the isotropic component U is positive for all values of q, and 
the lower bound on the anisotropy limit are s~ = 0.212675 (for 
7 = 1) and s~ = 0.0141 (for 7 = 2). These two limits are rep- 
resented as solid triangles in Fig. |2] As expected, their position is 
found between the NC and the WSC loci. Instead, in the 7 = 
case the function U is negative and the function V is positive for 
g < — 2.4: however sj > st , and the model is inconsistent. Note 
that the OM anisotropy limit for 7 models without dark halo, de- 
rived in C99 from their DF, is Sa (0) ~ 0.445, s~(l) ~ 0.202, 
and Sa(2) = 0, so we conclude that the presence of the DM halo 
slightly reduces the ability of the stellar density distribution to sus- 
tain radial anisotropy. 

In Fig.[3]the DFs of 7 = 1 and 7 = 2 models are presented, 
in the isotropic (solid), mildly anisotropic (dotted), and maximally 
anisotropic (dashed) cases. The dashed curves are very similar to 
the analogous curves in Ciotti & Lanzoni (1997, Fig. 2), and C99 
(Figs. 2 and 3), revealing the common qualitative behavior of OM 
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anisotropic DFs near the consistency limit, i.e. tlie fact that the in- 
consistency manifests itself in general at intermediate energies (see 
also Ciotti & Morganti 2008 for a discussion). 



4 JEANS EQUATIONS WITH OM ANISOTROPY 
4.1 Spatial velocity dispersion 

The solution of the spherical Jeans equations with general (radial or 
tangential) anisotropy has been given by Binney & Mamon (1982) 
and for OM systems is given by 



G 



2A( s) + sll{s) 
s^ + sl ' 



(29) 



where Sa = r^jr^. For the present models the explicit expression 
of the functions A and / are given in Appendix C for generic 7, 
however the resulting formulae are not particularly illuminating, 
as is common for this kind of models. Nonetheless, it is of some 
interest that in the 7 = case the velocity dispersion does not 
present "unreasonable" behaviour, even though we know that the 
model is physically inconsistent. Instead, the asymptotic analysis of 
(Jr at large radii and near the center provides helpful informations 
that will be used when discussing the projected velocity dispersion 
profile of the models. 

In the radial region r ^ r, , both A and J can be easily eva- 
luted, because the stellar density profile is asymptotic to the 
profile independently of 7. In any case tends to a constant: this 
is not surprising, as an elementary integration shows that the ve- 
locity dispersion profile of power-law densities in a total density 
profile r~'^ tends to a constant. In fact, an explicit calculation (or 
the expansion of eqs. [C1]-[C3]) shows that for s cxj 



2s' 



■4(s2 



(30) 



Therefore, in the fully isotropic case (Sa ^ 00) the radial veloc- 
ity dispersion at large radii is half of the model circular velocity. 
If some radial anisotropy is present, then at r ^ ra the orbital 
distribution becomes fully radially anisotropic, and accordingly the 
(square) intrinsic radial velocity dispersion increases by a factor of 
two when compared to the isotropic case. 

The situation is more delicate for r 0. In fact, from asymp- 
totic expansion of the integral in eq. l l29t it follows that the central 
behavior of (Tr is coincident with that of the isotropic case, and the 
product p,aj diverges for r ^ indipendently of the value of ra 
and 7. In addition, for ra > and 7 > 0, the product p.a^ di- 
verges as p,, so that converges to a finite value except for the 
7 = models: 



-logs, (7 = 0), 

-, (0<7<2). 

7 



(31) 



This is relevant from the modelistic point of view, as it is well 
known that self-gravitating isotropic 7 models present a depression 
of their velocity dispersion near the center with CTr(O) = (except 
for the 7 = and 7 = 2 models, see Bertin et al. 2002 for a general 
discussion of this phenomenon; see also Binney & Ossipkov 2001). 

Before discussing the projected velocity dispersion, it can be 
of interest in applications to have the analytical expression of the 
total kinetic energy of the stellar component. As is well known. 



from the virial theorem this quantity is independent of the spe- 
cific orbital anisotropy considered, and can be obtained without 
using the explicit solution of the Jeans equations. In fact 2if, = 
/ p*Tr((T2)d^x = / (x, V(^t)p* d^x = 47rG7^A/. A.rJ , where 
the last identity holds for any system of finite total mass M* in the 
gravitational field of the density distribution ([6]( (see also Kochanek 
1994), and in the present case 



GMl (3 - 7)7^ 



(32) 



Thus, if one defines the one-dimensional stellar virial velocity dis- 
persion as 3iV'/,(Tv/2 — Kt, it follows from the equation above 
that 



2 GM, (3 - 7)7^ vt 

f^v = 5 = -IT 

r, 3 3 



independently of the value of 7. 



(33) 



4.2 Stability 

Equations (29) and (32) can be used to obtain indications about 
the minimum admissible value of Sa as a function of 7 to pre- 
vent the onset of radial orbit instabihty. A complete stabihty 
analysis is beyond the task of this work, requiring N-body sim- 
ulations or normal mode analysis, but some interesting con- 
clusions can be equally derived from the work of Fridman & 
Polyachenko (1984). These authors argued that a quantitative 
indication on the maximum amount of radial anisotropy sus- 
tainable by a specific density profile is given by the stability 
parameter ^ = 2i^r*/A't,, where A'r, and Kt„ — K, — A'r, are 
the radial and tangential component of the kinetic energy ten- 
sor, respectively. From its deflnition ^ ^ 1 for Sa — > 00 (glob- 
ally isotropic models), while ^ ^ 00 for Sa — > (fully radially 
anisotropic models), and for one-component systems stability 
is associated with the empirical requirement that ^ < = 
1.7 ± 0.25; the exact value of is model dependent (see, e.g., 
Merritt & Aguilar 1985; Bertin & Stiavelli 1989; Saha 1991, 
1992; Bertin et al. 1994; Meza & Zamorano 1997; Nipoti, Lon- 
drillo & Ciotti 2002). Here we are considering two-component 
systems, however N-body simulations have shown that the pres- 
ence of a halo does not change very much the situation with re- 
spect to the one-component systems (e.g., see Stiavelli & Sparke 
1991, Nipoti et al. 2002). Therefore, in the following discussion 
we assume as a fiducial maximum value for stability = 1-7. 

The parameter ^ for the present models is independent of 
TL, and it cannot be expressed by using elementary functions, 
so that we explore its value numerically. In Fig.|4]we plot ^ as 
a function of Sa for 7 = 1 /2, 1 and 2, and the asymptotic flat- 
tening to unity for increasing isotropy is evident. It is appar- 
ent that the stability criterion requires minimum anisotropy 
radii appreciably larger than those obtained from the consis- 
tency analysis (see Sect. 3.2). In addition, stable stellar distri- 
butions with shallower central density profile require more and 
more isotropic velocity dispersion, confirming the trend already 
found for one and two-component 7-models (e.g., Carollo et al. 
1995; Ciotti 1996, 1999). So, it is likely that the more radially 
anisotropic models with positive DF are prone to radial orbit 
instability. 
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Figure 4. The stability indicator § as a function of the normalized 
anisotropy radius Sa = fa/r,, for models with 7 = 1/2, 1, 2 (dashed, 
dotted, and solid lines, respectively). The horizontal line marks the fiducial 
stability limits = 1-7. 



4.3 Projected velocity dispersion 

The projected velocity dispersion associated with a general 
anisotropy function /3(r) (see eq. 1171 ) is given by 



E.(i?)a^(i?) = 2 



r,2 



p*{r)a^{r)r 
Vr2 - i?2 



dr. 



(34) 



Unfortunately the projection integral above cannot be evaluated an- 
alytically for generic 7 in terms of elementary functions. However, 
as for the projected stellar density, interesting informations can be 
obtained in two relevant radial regions, i.e., outside the core radius 
and near the center. In practice, the external regions are defined as 
the radial interval where the stellar density profile can be approx- 
imated as a pure power-law of slope —4. In this region the pro- 
jection integral can be evaluated for generic values of Sa and the 
asymptotic result is 



aliR) 



i/2 . 



4si( 



si + f) 



,2-13/2 



(35) 



where 77 = i?/r*. In the fully isotropic case the dimensionless ra- 
tio (Tp [R) /vc tends to 1/2, so that the projected velocity dispersion 
coincides with the (constant) isotropic velocity dispersion (as ex- 
pected), while in the completely radially anisotropic case (or for 
f) 3> Sa) the ratio converges to 1/ \/8- 

The case of the central regions is more complicated. In fact, 
integral l |34| l for R ^ Q converges when < 7 < 1, and diverges 
for 1 < 7 < 2, as can be easily proved by using eqs. (C5)-(C6). In 
the divergent case both the projection integral and the projected sur- 
face density E, are asymptotically dominated by their integrands 
for r (as is intuitive, the cusp dominates and the contribution 



of foreground stars and background stars is negligible), and o-p(R) 
can be properly defined as the limit value of their ratio for R 0. 
A simple calculation shows that 



ap(0) = a,(0) 



1 < 7 < 2. 



(36) 



In other words, for the models with the centrally divergent sur- 
face (stellar) density, the projected central velocity dispersion co- 
incides with the central radial component of the isotropic velocity 
dispersion. Thus it does not depend on the model anisotropy radius, 
at variance with the projected velocity dispersion in the external 
galactic regions. 

In the convergent case < 7 < 1 the central value of the 
projection integral depends on the whole profile of the integrand, 
therefore the projected central velocity dispersion depends also on 
Sa. For generic 7 and Sa, the quantity crp(O) is expressible in terms 
of hypergeometric 2F\ functions. However a simple form of the 
projection integral evaluated at i? = 0, useful in numerical inte- 
grations, can be obtained by inverting the order of integration in 
eq. (O: 



(1 + ^)*- 



Sa \ S , 

1 H — ^ I arctan — ds. 

Sa 



(37) 



In the fully isotropic case (Sa 00) the projection integral can be 
evaluated analytically and the result is 



ap(0) 



< 7 < 1, 



(38) 



which is independent of 7. This is not surprising, as it is easy to 
show, by inverting order of integration, that for any density profile 
with finite central projected density, isotropic orbital distribution, 
in the total gravitational field produced by density distribution (6), 
identity (38) holds. 

Equations l l35t , l l36t . and l l38t open the interesting possibility 
to consider the ratio of the outer to the central projected velocity 
dispersion, a quantity that can be expressed in a very simple way. 
In particular, the ratio depends on the shape of the stellar density 
slope 7, on the outer observational point R, and finally on Sa. Thus, 
at least in principle, for galaxies well described by a 7-model im- 
mersed in a total density profile oc r~^, it could be possible to de- 
termine Sa from observations, assuming OM anisotropy. In Fig. 4 
the ratio is plotted as a function of the anisotropy radius for the 
representative values of the density slope 2, 1, and 1/2 (in this lat- 
ter case the WSC limit for consistency is Sa>0.45). All the ex- 
pected trends are apparent, iin particular the decrease of the ratio 
(jp(i?) /crp(O) for decreasing Sa. This is due, for 1 < 7 < 2, to the 
decrease of o-p{R) in the external regions of radially anisotropic 
models. In the 7 = 1/2 case, reported as an example of models 
with central slope in the range < 7 < 1, the larger decrease is 
due to the additional effect of the increase of o"p(0) for decreasing 
Sa. Overall, the kinematical ratio Op{R) /crp(O) for mildly-strongly 
anisotropic models (i.e., Sa<l) is ~ 15% — 20% lower than in the 
corresponding isotropic cases. 

4.4 Some additional considerations 

For sake of completeness, we summarize some additional results 
on velocity dispersion. For example, the central velocity dispersion 
of 7 models in the presence of a black hole can be found in the 
literature. Here we just recall that the velocity dispersion diverges 
for r ^ as r~^'^^ (e.g., see C96, Baes & Dejonghe 2004, Baes et 
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al. 2005; for the case of oblate 7 models, or two-component oblate 
power-law models with central black hole see also Riciputi et al. 
2005, Ciotti & Bertin 2005). It follows that in the present context 
a central black hole would produce an identical kinematical signa- 
ture, as sufficiently near the center the total mass is fully dominated 
by the black hole. 

As a second case, we consider the spatial and projected veloc- 
ity dispersion of a two-component galaxy model made by the su- 
perposition of a stellar component described by a 7 model (where 
for simplicity we restrict to the interval 1 < 7 < 2), and a dark 
halo component described by a Jaffe model (7 = 2). The interest 
of these models is due to the fact that in the inner regions they be- 
have as the models subject of this work, but in the external regions a 
Keplerian decline is present, due to the halo finite total mass. From 
eq. ([TJ, the Jaffe profile of total mass A/h = TZM^ and scale length 
rii — l3r, can be written as 



Ph(r) 



(3-7) s2(/3 + s)2 



(39) 



The Jeans equations for this class of models cannot be solved ex- 
plicitly in terms of elementary functions for generic 7, even though 
special explicit cases can be easily found (e.g., see Ciotti et al. 
1996). For this reason we just evaluate the asymptotic leading term 
of the velocity dispersion for s — > 00, obtaining 



2 GM,{1 + TZ) 5s^+3sl 
^ r, 15s(s2 + s2)' 

and its projection through eq. l l34t . giving 

GM, 4(1 +7^) 



(40) 



ISttt; 



2 + 



^2)3/2 



(41) 



where 77 = R/rt - The expression in square parentheses converges 
to 2 in the isotropic case. 

The behavior of the velocity dispersion in the central regions 
requires some additional discussion. In fact, for 1 < 7 < 2 the 
central projected velocity dispersion will depend on the central, 
isotropic spatial velocity dispersion only (for the reasons described 
in Section 4.2, and excluding the purely radial case). In addition, 
the self-contribution to the central velocity dispersion of the 7 
model is zero (except for the 7 = 2 case, e.g., see Bertin et al. 
2002). Summarizing, the projected central velocity dispersion of 
these models is coincident with the isotropic spatial velocity dis- 
persion of 7 models in the presence of the Jaffe halo only (with the 
exception of the 7 = 2 case). Therefore, 



a^(0) = a,2(0) = 



CM. 



2(3 ' 



(7 = 2). 



(42) 



4.5 Velocity profile 

Not surprisingly, the velocity profile VP cannot be expressed in 
terms of elementary functions, however acceptably simple formu- 
lae can be obtained at large radii and in the central galactic regions. 
The starting point is to consider the normalization of eq. ( I20t . given 
by 



E.VP : 



g{s,r])sG{qm)ds 



(43) 




Figure 5. The ratio ap{R)/crp{0) as a function of the normalized 
anisotropy radius Sa = ra/r*, for models with 7 = 1/2, 1, and 2. Sohd 
lines refer to R = 2R.c. dotted lines to R = ARc- 



where q„ 

G{qn.) = 



dg 



and 



(44) 



At the very center of the stellar system, where the stellar density 
profile is proportional to s^^, it is not difficult to show that for 

1 < 7 < 2 



VP . 



27rt)2 



(45) 



independently of the (positive) value of the anisotropy radius. In- 
stead, for < 7 < 1, as for the projected velocity dispersion a 
numerical integration is required. 

In the external galactic regions, where the stellar density pro- 
file can be approximated as power-law, the VP can be written 
as a quite simple integral, depending on the dimensionless ratios 
v\Jvc and r;/sa, that can be evaluated numerically. Here we just 
report the formula for the isotropic case, where 



(46) 



In both the reported formulae, the Gaussian signature is apparent. 



5 CONCLUSIONS 

In this paper a family of spherical, two-component galaxy models 
with a stellar density profile described by a 7 model and a total 
(stars plus dark matter) density profile oc r^^ at all radii has been 
investigated, under the assumption that the internal dynamics of the 
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stellar component is described by Osipkov-Merritt anisotropy. The 
models are fully determined when the inner density slope 7 and the 
anisotropy radius ra of the stellar component are assigned, together 
with a density scale for the total density profile. The dark matter 
halo remains defined as the difference between the total and stellar 
density profiles. The main results can be summarized as follows. 

• After having provided the most common structural quantities 
of the models that are of interest for observations, limitations on the 
total density scale as a function of 7 are analytically determined by 
requiring the positivity and monotonicity of the dark matter halo 
distribution. In particular, the request of positivity limits the range 
of acceptable stellar density slopes to < 7 < 2. Models corre- 
sponding to the minimum total density scale (for given 7) are called 
minimum halo models. The central density profile of the dark mat- 
ter halo diverges as in general, but in the minimum halo 7 = 2 
model (in which the positivity and monotonicity limits coincide), 
the central dark matter profile is oc r^^. 

• The minimum value of anisotropy radius corresponding to a 
dynamically consistent stellar component has been derived analyt- 
ically as a function of 7 by using the necessary and sufficient con- 
ditions of CP92. As expected, an increase of 7 results in a decrease 
of the minimum value of the anisotropy radius ra required by con- 
sistency. It is also proved that models with 7 > 4/17 are certainly 
consistent for sufficiently isotropic velocity dispersion, while for 
centrally shallower stellar density profiles the sufficient condition 
for consistency is never satisfied. The necessary and sufficient con- 
dition for the halo consistency are also analytically obtained, and 
the minimum halo models corresponding to (isotropic) dark mat- 
ter halos are derived. In the case 7 = 2, the minimum halo co- 
incides with the minimum halo obtained from the positivity and 
monotonicity conditions. 

• The phase-space DF of the stellar component for 7 — 0, 1, 2 
is analytically recovered in terms of Polylogarithms and exponen- 
tials. It is found that the 7 = model is inconsistent no matter 
how much anisotropy is considered. Instead, the isotropic 7 = 1,2 
models have a positive DF, and the true critical anisotropy ra- 
dius for consistency can be determined directly from their DF. A 
comparison with the analogous study of one-component 7 models 
shows that the presence of the halo sligthly reduces the maximum 
amount of sustainable radial anisotropy. The obtained values of the 
anisotropy radius are independent of the total density scale TZ. 

• The Jeans equations for the stellar component are solved ex- 
plicitely for generic values of 7 and ra in terms of elementary func- 
tions. The asymptotic expansions of ar for r ^ and r ^ 00 
are obtained, and it is shown that ar tends to finite (non-zero) val- 
ues (except for the divergent central velocity dispersion of 7 = 
model) which are simply related to the model circular velocity. The 
projected velocity dispersion <Jp{0) cannot be calculated analyti- 
cally, in general. However, by asymptotic expansion of the pro- 
jection integral, exact values at large radii and at the center are 
obtained. In particular, it is shown that for 7 > 1, and indepen- 
dently of the value of the anisotropy radius, crp(O) coincides with 
the central velocity dispersion (7r(0) in the isotropic case. Instead, 
for < 7 < 1, (Tp(0) depends also on Sa. In the anisotropic 
case crp(O) cannot be obtained analytically; however, a very simple 
form of the projection integral suitable for numerical integrations 
is given. In the isotropic case the integral can be evaluated analyt- 
ically, and o-p(O), as expected, coincides with the model circular 
velocity. 

• Finally, we have shown that the Velocity Profile of the models 
can be obtained in a very simple form (Gaussian) near the center 



(independently of the value of the anisotropy radius and for 1 < 
7 < 2), and at large radii (in the isotropic case). 

We conclude by noting that these models, albeit highly ide- 
alized, seem to suggest two interesting remarks of observational 
character. The first is that in real galaxies with a total r^^ density 
profile and sufficiently peaked stellar density (i.e. 7 > 1), mea- 
sures of central velocity dispersion should not be strongly affected 
by radial anisotropy. Second, for given central stellar density slope 
7, measures of the projected velocity dispersion at the center and in 
the external regions are able, at least in principle, to determine the 
value of the anisotropy radius under the assumption of Osipkov- 
Merritt anisotropy. 



ACKNOWLEDGMENTS 

L.C. wishes to thank Carlo Nipoti and Leon Koopmans for use- 
ful comments. We also thank an anonymous referee for useful 
comments. 



APPENDIX A: NECESSARY AND SUFFICIENT 
CONDITIONS FOR MODEL CONSISTENCY 

The condition for the positivity of the halo density profile ph can 
be easily established from eq. dlOb . In fact, for < 7 < 2, 7?. must 
be greater or equal to the maximum of the radial function inside 
the parentheses, and simple algebra shows that the maximum is 
attained for 
2-7 



(Al) 

in particular, for the Jaffe model the critical point is reached at the 
center. The monotonicity condition is obtained requiring that the 
radial derivative of ph is nowhere positive, and this happens if and 
only if 



7^ > 



8^-^(7 + 4s) 



2(1 



Vs > 0. 



(A2) 



Thus, TZ must be greater than or equal to the maximum of the radial 
function above, that is reached at 



12 - 77 + ^(4 -7) (36- 177) 



16 



(A3) 



and again for 7 = 2 the maximum is reached at the center. Finally, 
the application of the WSC to an isotropic halo in order to have 
phase-space consistency is given by eq. i24i with ra 00, so that 
g = ph- The condition becomes 



TZ > 



[16s^ + (97 - 4)s + 7^] 
4(1 + 5)*^-^ 



Vs > 0. 



(A4) 



The study of the maximum of the rh.s. of equation above leads 
to a cubic equation. We do not report here the solution Sm(7) 
corresponding to the maximum, as it can be easily obtained, but 
for reference we just report three special values Sm(0) ~ 2.2049 
Sm(l) — 1.2079 and Sm(2) = 0, corresponding to TZm{0) ~ 
0.07735, 7^m(l) ~ 0.17487, and TZm{2) = 1, respectively. 

The application of the WSC to the stellar component reduces 
instead to the study of 



16s^ + (97-4)s+7 +s' 



4s2 + (57-12)s + (7-2) = 



> 0, (A5) 
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so that the function F in eq. l l25t is a quadratic polynomial, and the 
determination of the sets C+ and C_ is straigthforward. In partic- 
ular, C- is not empty for < 7 < 4/17, and the function G is 
positive there. 



APPENDIX B: PHASE-SPACE DF FOR THE STELLAR 
COMPONENT 

Here we give the explicit evaluation of the DF for the three integer 
values 7 = 0, 1, 2. In fact, in these cases one can change the vari- 
able integration by defining t = ^/W~~q, so that e* = e' and 
the integration interval is mapped into (0, 00). Expansion in simple 
fractions, factorization of outside the integrals and repeated dif- 
ferentiation with respect to e"'' under the sign of the integral finally 
shows that for 7 = 

[Li-9/2(y) - 6Li_7/2(y) -I- llLi__5/2(y) - 6Li_3/2(y)] , (Bl) 

V{q) = ^ X [Li-9/2(y) - Li_5/2(y)] , (B2) 

where y = —e~'', Lis{z) — z^{z, s, 1) is the so-called Polyloga- 
rithm function, $(2, s,a) is the Lerch function and dLis (z) /dz — 
Lis-i{z)/ z (e.g., Erdelyi et al. 1953). A similar treatement for the 
7=1 case gives 

+ U_r/2{y) ~ 5Li_5/2(y) + 6Li_3/2(y)] , (B3) 



V{q) = ^ [Li-7/2(y) - Li_5/2(y)] 



(B4) 



and finally for 7 = 2 

U{q) = V^x [2^/^e~^^ ~ 2e-« + Li_5/2(y) - 3Li_3/2(y)] , (B5) 
V{q) = X [Li_5/2(y) - Li_3/2(y)] . (B6) 

APPENDIX C: VELOCITY DISPERSIONS 

The isotropic function / in eq. ( |29l > for 7 7^ 0, 1, 2 is given by 

r_6s^ + 6(3 - 7)s2 + 3(3 - 7)(2 - 7)^ + (3 - 7)(2 - 7)(1 " 7) 



57(1 + s)3-7(3 - 7)(2 - 7)(1 - 7)7 
6 



4s4' 2s2' 
while for s — * 

log s, 7 = 



(3-7)(2-7)(l-7)7' 



while 

1 + s ll + 15s + 6s2 
s 6(1 + s)3 

1 5 -I- 4s o, 1 + s , 

< ~ + ^77TT~^ " 3 log , 7 = 1 

s 2(1 + s)^ s 

1 + s , l-3s~6s^ 

2.^(1 + .) ' -^^2. 

The function A is given by 

1 (3-7 + s)s^~^ 

(3 - 7)(2 - 7) (3 - 7)(2 - 7)(1 + s)^-^ ' 



(C2) 



log 



1 + s 



(C3) 



s 1 + s ' 
For s 00 and < 7 < 2 



7 = 2. 




and 



A~<^ (3_^)(2_7) 
I - log s, 7 = 2 



, < 7 < 2 



(C4) 



(C5) 



(C6) 



(CI) 
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